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S3 

d [ Abstract. The pseudo-spin symmetry is reviewed. A mapping that 

produces the separation of the total angular momentum into pseudo- 
orbital and pseudo-spin degrees of freedom is discussed, together with 
the analytic transformations that take us from the normal parity space 
to the eigenstates of a pseudo-oscillator with one quanta less. The many- 
particle version of the unitary transformation to the pseudo-SU(3) space 
I is established. As an example, these symmetries are used to describe the 

double beta decay phenomenon in heavy deformed nuclei. 



o 

INTRODUCTION 

The pseudo-spin symmetry was introduced in [1,2], where a separation 
"q ■ of the single particle Hilbert space into pseudo orbital and pseudo-spin de- 

grees of freedom was proposed. This new coupling scheme provides a sim- 
ple interpretation of many striking features of the Nilsson orbits: i) The 
observed approximate degeneracy of the single particle levels of the type 
(N, I, s)j and (N, I + 2, s)j + 1 for the spherical case or the one associated 
to the levels with asymptotic quantum numbers [N,N Z ,A] Q = A + 1/2 and 
[N, N z , A + 2}Q = A + 3/2 for the large deformation case. These energy orbitals 
can be considered as pseudo-spin orbit doublets, implying that the strength 
of the pseudo spin-orbit interaction is small, ii) The expectation value of the 
pseudo-spin operator, iii) the calculation of the decoupling parameters, and iv) 
the matrix elements of the Coriolis interaction [1]. More recently the pseudo- 
spin symmetry has been considered to be useful to characterize identical bands 
[3] and as a signature of super deformation [4]. 

For heavy nuclei the nuclear shells can be divided into two parts: the asso- 
ciated to the single particle orbits of the same parity, which are called normal 
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parity levels and the intruder that comes from the shell above of opposite par- 
ity and thus it has been called unique or abnormal parity orbital. Thus for a 
given shell, N a one has j£ = {1/2, 3/2, • • • , N a - 1/2 }, j£ = N a + 3/2, where 
a = n or v for protons and neutrons, respectively. This separation of the space 
for the single particle states of Nilsson or Wood Saxon potentials is reasonable 
because it has been proved that the intruder orbitals are nearly unmixed for 
standard deformations [5]. When relabelled in the pseudo spin scheme, the 
normal parity levels form a major shell for a pseudo-oscillator potential with 
N a = N a — 1 quanta. Since the symmetry of this oscillator is of course SU(3), 
one may consider the pseudo-SU(3) coupling scheme. 

In the Second Section, we review the transformation to the pseudo space 
introduced by Bohr, Mottelson and Hamamoto (BHM) [6] and the mappings 
from oscillators to pseudo-oscillators wave functions for the spherical and 
asymptotic cases [7]. We emphasize the difference between the pseudo-spin 
and the pseudo-oscillator symmetries. The first one implies the separation 
of the normal parity degrees of freedom into pseudo spin and pseudo orbital 
spaces (called pseudo-space) while the second one for the one particle case 
is related with the use of oscillator wavefunctions with spin to describe the 
pseudo-space. At the end of this section, the many-particle unitary trans- 
formation to the pseudo-SU(3) space is introduced and applied to one body 
operators. In the final section, theoretical calculations for the 2v double beta 
decays of 150 Nd and 238 U are presented [8]. 

UNITARY TRANSFORMATIONS 

The BHM unitary transformation acts on the angular-spin parts of the 
wave functions of the particle, without affecting the radial motion. It is given 

— # 

by the scalar product between the spin operator S and f a unit vector in the 
direction of the position of the particle, Ubhm — 2ir ■ S. It is straightfor- 
ward, to calculate the action on the single-particle wave function of the three 
dimensional harmonic oscillator with quantum numbers 

1pN(l,s)jm(r) = Rm(r) J2 SCr \j m ) Y l^: ^Xsa , (1) 

which can be written as a product of the radial part times the superposition 
of two spinor wave functions [6]. This transformation gives rise to the results 
I = j ±1/2 i— > I — j =p 1/2 . Thus it produces the breaking of the normal 
degrees of freedom into pseudo-spin and pseudo-orbital spaces, but it does not 
give the wave functions of a pseudo-oscillator with one quanta less. 

In the set of wave functions (2) corresponding to a given N value, we disre- 
gard those with j = N+^, which define the unique or abnormal space, and only 
consider the remaining wave functions of the normal parity orbitals. These 



states can be mapped onto the eigenstates \N(ls)jm) of a pseudo-oscillator 
via the following unitary operator [7] 

U = 2(g.S)(N-2L.S)- 1 *, (2) 

— * 

where £ is the annihilation harmonic oscillator quantum. The action of U onto 
eigenstates of a three dimensional harmonic oscillator gives the result 

U\N(l,s)jm) = \N(ls)jm), (3) 

with the relations between the labels N — N — 1, s — s, j — j, m — to, and 
I — I ± 1 according to whether j = I ±1/2 . Although the previous mapping 
connecting the normal parity harmonic oscillator eigenstates with the full set 
of eigenstates of a pseudo-oscillator was proposed long ago [1, 2], the explicit 
form of the unitary operator (3) was constructed for the first time in [7] . More 
recently, the transformation to the pseudo-oscillator space has been carried 
out by using the Algebraic Generator Coordinate Method [9]. 

Now we consider the asymptotic wave functions of the Nilsson Hamiltonian 
[7] associated with large deformations. Then the Nilsson orbitals are described 
by the states characterized by | [N, N z , A]Q), the cylindrical harmonic oscil- 
lator states with spin. The degeneracy mentioned above can be explained 
through the unitary transformation 

C/oo = 2 (£ + S + + £_S_) (N p - 2L Z • , (4) 

where the operators £± and r]± are the spherical components of the creation 
fj and annihilation £ harmonic oscillator operators, the spin operators are 
expressed in terms of the Pauli matrices S m = a m /2 and N p is the number 
operator in the plane. 

The Uoo and their corresponding hermitean conjugated operator are 
unitary, if they act onto states belonging to the normal space orbitals, that is 
if we disregard the levels of a given iV shell that satisfy = N p + 1/2. Thus, 
the unitary transformations for the spherical and asymptotic Nilsson orbitals 
to the pseudo-oscillator wave functions are not equivalent. This fact which 
has not been previously emphasized, is reflected in the separation itself into 
normal and abnormal parity spaces. However for standard deformations the 
exact wave functions of the Nilsson Hamiltonian can be evaluated and it is 
observed that there is no mixing of the unique (spherical) parity orbital with 
the remaining levels of the shell [10], which suggest that for these deformations 
the spherical transformation is the most appropriate. The division of the shell 
model space into normal and unique as proposed in the asymptotic limit may 
be useful in applications to superdeformation phenomena. 

For a system of n particles the unitary transformation to the pseudo- 
oscillator space is defined by 

n 

U = UU.. (5) 



In the Fock space, if we denote by U the corresponding unitary operator 
associated to U, then a general one body operator F is mapped to the pseudo- 
SU(3) space by F = UFU^ . If two single particle state vectors |a) and \a) are 
related by a unitary transformation U (see Eq. |2|), then the fermion operators 
are related by Uaj^W = a' & . Then, F in the pseudo-space is given by 

F= E ' E ' (N^^hf^FxlN^hj^a^,^ a^.-, (6) 

N,l,j,/j, N',l',j',n' 

where the primes on the sum indicate that the unique parity orbitals are 
excluded. A similar result is obtained for the two body operators, that is only 
the labels of the fermion creation and annihilation operators are changed by 
its pseudo quantum numbers. The previous results give formal support to the 
recipe indicated in Refs. [1, 11]. The many body expressions in the pseudo 
space of the SU (3) generators are given by the series expansions 

L g = k L L q + • • • , g M = k Q Q M + • • • , (7) 

where only the leading terms of the series are indicated. These operators Lq.Q^ 
have the same form as the 577(3) generators but of a shell with one quanta 
less, but we want to emphasize that they are not the transformed operators 
to the pseudo space. 

DOUBLE BETA DECAY 

Double beta decay is a rare transition between two nuclei with the same 
mass number A involving change of the nuclear charge number Z by two 
units. This exotic phenomenon is a useful tool to test the lepton number 
conservation, neutrino properties and models of nuclear structure [12] . It can 
be classified into various modes according to the light particles besides the 
electrons associated with the decay. The two neutrino mode (ftfav), in which 
two electrons and two neutrinos are emitted, takes place independently of the 
neutrino properties, and conserves the electric charge and lepton number. The 
Ou mode violates lepton number conservation and therefore it is forbidden in 
the standard electroweak theory. To proceed the Ou decay, the virtual neutrino 
must be emitted in one vertex and absorbed in the other one, thus it is required 
that: i) the exchanged neutrino is a Majorana particle {y = v) and ii) both 
neutrinos have a common helicity component [13]. 

Next, we restrict to describe, within the pseudo-SU(3) formalism, the 2v 
double beta decays of 150 Nd and 238 U. The decay rate of the 2i/-mode can be 
calculated through the formulae 

(rli 2 )- 1 = G 2v \ Mg T \\ (8) 



where G 2 v is a kinematic factor and M^J = M 2v is a nuclear matrix element 
strongly dependent on the considered nuclear model, that is 

m 2 , = Ett— ^ — W V$ l|r|| 1+) (i+ ||r|| o+>. (9) 

The T is denoting the Gamow- Teller operator T m = J2 S Cms^J and the E = 
\ Q/3/3 + n2 e c 2 is the half of the total energy released. The E N gives the energy 
of the intermediate state | ljv) - The E t is the energy of the ground state of 
the initial nucleus |0+) and the ket |0j) describes the ground state of the final 
nucleus. 

In order to compute M 2u , we have to perform a sum over all the intermedi- 
ate states . Fortunately, an alternative form of calculate this matrix element 
has been developed [14], i.e., 

M ^ = 4-<o? I £ (-i) m r_ m F m |o+), (io) 

where the operator F m is defined by: 

oo MU 

F m = Y. -TT- [H,..., [H, T m ] . . .]<*-**»«> . (11) 

A reasonable model for describing spectra and BE2 transitions of heavy de- 
formed nuclei is [8]: 

H = E H a - \ X Q a ' Q a + Ci K 2 + C 2 L 2 , (12) 

where H a denotes the spherical Nilsson Hamiltonian for neutrons or protons 
plus a constant term V a , which represent the depth of the potential well. 
The quadrupole-quadrupole interaction in a given shell and K 2 is a linear 
combination of L 2 , X3 and X4, which are rotational scalar operators built 
with generators of the algebra of SU(3) [15]. 

Notice that the quadrupole-quadrupole force, L 2 and the K 2 interaction 
are independent of the spin degrees of freedom and symmetric in the neutron 
and proton components. To evaluate (13), we express H a and T m in the second 
quantization formalism, 

H a = hw £ e a (77, /, j) aj, 1 jma a vl i jma , 

rjljm 

T m = <r(7T, I/)A(7T, 1/, m) , 

7T V 

with A(ir, v, m) = [a' l 1 iyJL denoting the angular coupling of pro- 
ton creation and neutron annihilation operators, e a the single particle energies 



(13) 
(14) 



and 

<r{*, v ) = H^-^^- (VnK ^'JttIWWv Ju^;ju) ■ (15) 
Afterwards some algebraic manipulations, M 2v is rewritten by 

jij V- V- / / ,\ ( \ < / | A(n',u')-A(n,i<) | 0+) 

M 2 ,= £„„£„,,,, <t(tt ,1/ ) a(ir,v) Eq f + ) _ ) (16) 

In the pseudo SU(3) Hilberty space the above sum is restricted to r) n = r) v , 
K = lui Vn' = f]v'iK' — Iv'i j-n — ju + 1- Following [8], the single particle 
energy difference in the denominator takes the form 

eivJJn) -z(v],l,3v) = -huk w 2j n + A c ■ (17) 

The constants k a are well known [16] and Ac is used to determine the value 
V v — V n . The Ac is the difference Coulomb energy and it is evaluated by the 
expression 

A c = [2Z + 1 - 0.76((Z + I) 4 / 3 - Z A ^)]MeV . (18) 

The description of the correlated deformed ground states is done using the 
pseudo SU(3) scheme for the normal parity space and seniority zero configura- 
tions, that is all the nucleons coupled by pairs to angular momentum zero, for 
the unique part [11, 15]. The occupancies of these spaces are determined from 
the corresponding Nilsson diagrams by selecting a reasonable deformation and 
filling each level with a pair of particles in order of increasing energy. These 
numbers fix the totally antisymmetric irreducible representations (irreps) of 
the unitary groups associated to the normal U((N a + l)(N a + 2)) and unique 
U(2N a + 4) spaces. 

For the normal parity space, one separates the degrees of freedom in 
pseudo-orbital U(^) and pseudo-spin U a (2) parts, with Q% = (N a + l)(N a + 
2)/2 and their irreps {f a } are characterized by the partitions of the num- 
ber of particles in the normal part, ■ Then of the pseudo SU(3) irreps, 
(A Q ,/i Q ), contained in {f a }, one considers those with maximum eigenvalue of 
the Casimir operator, (C 2 ) a = (X a + /i Q + 3) (X a + ji a ) — \ a (J, a . In Table 1, 
the results found for the participant nuclei in the double beta decays of 150 Nd 
and 238 U are presented. 

Table 1. Ground states in the Pseudo-SU(3) coupling scheme. 



NUCLEUS 


< 










U(S1») 


SU n (3) 


5^(3) 


150 Nd 


6 


4 


6 


2 


{2 3 } 


{23} 


(12, 0) 


(18, 0) 


150 Sm 


6 


6 


4 


2 


{23} 


m 


(12, 0) 


(12, 2) 


238 jj 


6 


4 


12 


8 


{23} 


{2 6 } 


(18, 0) 


(36, 0) 


238 p u 


6 


6 


10 


8 


{23} 


m 


(18, 0) 


(30, 4) 



Table 2. Theoretical estimates for the half-life /3/3-decay in the 1v mode. 



Transition 


<o+|r 2 |o: 


h > E[MeV] 




OH13] 


l50 Nd _^150 Sm 


1.31 


12.2 


6.0 x 10 18 


9-17 x 10 18 


238 jj _238 p u 


1.51 


16.8 


1.4 x 10 21 


2 x 10 21 



Finally we used the strong coupled limit [11], and from the Kronecker 
product (\ n ,/J, n ) x (A„,//„), the (\ n + \ u , /i n + fi u ) irrep will dominate the 
low-lying energy structure. 

In Table 2, we display the calculated Gamow- Teller matrix elements, energy 
denominators, predicted and experimentally determined double beta half lives 
of 150 Nd and 238 U. They are given in the last column of Table 2 and these half 
lives are lower bounds because the nuclear matrix elements could be reduced by 
considering, for example the inclusion of the pairing interaction, and therefore 
giving longer (never shorter) (3(3 half lives. 
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